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Abstract 


The problem of electromagnetic (EM) scattering from an aperture backed by 
a rectangular cavity recessed in a 3D conducting body is analyzed using the coupled field 
integral equation approach. Using the free-space Green’s function, EM fields scattered out- 
side the cavity are determined in terms of 1) an equivalent electric surface current density 
flowing on the 3D conducting surface of the object including the cavity aperture and 2) an 
equivalent magnetic surface current density flowing over the aperture only. The EM fields 
inside the cavity are determined using the waveguide modal expansion functions. Making 
the total tangential electric and magnetic fields across the aperture continuous and subjecting 
the total tangential electric field on the outer conducting 3D surface of the object to zero, a set 
of coupled integral equations is obtained. The equivalent electric and magnetic surface cur- 
rents are then obtained by solving the coupled integral equation using the Method of 
Moments (MoM). The numerical results on scattering from rectangular cavities embedded 
in various 3-D objects are compared with the results obtained by other numerical techniques. 

1. INTRODUCTION 

Electromagnetic scattering characteristic of metallic cavities is useful in studying 
radar cross section and electromagnetic penetration properties of objects consisting of these 
cavities as substructures. A large amount of analytical work has been done to characterize these 
cavity structures. A few references, but not a complete list, are given in [1- 8]. However, in 
these and similar work, it is assumed that the aperture backed by a cavity is in an infinite flat 
ground plane. For a characterization of an aperture formed by a cavity recessed in a finite gound 
plane or no ground plane, the asymptotic techniques described in [9-11] may be used. However, 


6 



these asymptotic techniques are applicable when the frequency is high or when the cavities are 
large in size compared to the operating wavelength. For cavities with size comparable with a 
wavelength, a rigorous integral equation formulation has been used to analyze cylindrical circular 
cavities [12-13]. However, the approach described in [12-13] uses entire domain expansion 
functions to represent surface current density and hence is limited to only cylindrical circular 
cylinders, A need, therefore, exists to develop general analytical tools to determine low 
frequency electromagnetic characteristics of open-ended waveguide cavities without a flang or 
recessed in 3D counducting surface. 

In this paper the problem of EM scattering of plane waves by a cylindrical cavity 
recessed in a 3D metallic object is studied. Using the equivalence principle, the electromagnetic 
field scattered outside the object is determined using free space Green’s function and the 
equivalent electric and magnetic surface currents assumed to be present on the outer surface of the 
object. The equivalent electric surface current is assumed to be flowing over the complete 3D 
surface including the aperture and the equivalent magnetic surface current is assumed to be 
flowing over only aperture. The field inside the cavity are obtained using waveguide modal 
expansion functions. Making the total tangential electric and magnetic fields across the aperture 
continuous and the total tangential electric field zero over the conducting surface only, a set of 
coupled integral equations is obtained. Expanding the surface currents in triangular subdomain 
functions [14] and using the Method of Moments, the coupled integral equations are reduced to 
algebric equations which are solved for the surface current densities. From the surface currents 
the radar cross sections of these cavities recessed in an arbitrarily shaped conducting objects are 
determined. For future reference, this method is refered to as Modal Expansion and Method of 
Moments ( ME-MoM). 
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The remainder of report is organized as follows. The formulation of the problem in 
terms of coupled integral equations using surface equivalence principle is developed in section 2. 
Numerical results on radar cross section of open-ended rectangular cavities recessed in a 
rectangular and finite circular cylinder are presented in section 3. Comparison of the numerical 
results obtained by the present method with other numerical techniques is also presented in 
section 3. The advantages and limitations of the present formulation are discussed in section 4. 

2. THEORY 


Consider a time harmonic electromagnetic plane wave incident on a 3D conducting 
object with an aperture backed by a rectangular cavity as shown in figure 1. The cavity is formed 
by a shorting plate at z = -L. To facilitate the solution of the problem, the equivalence principle is 
applied by using the equivalent surface currents as shown in figure 2. For determining the fields 

outside the cavity (exterior problem), we consider the equivalent currents ) and m a radiating in 

free space. The electromagnetic field inside the cavity ( interior problem ) is obtained using 
modal expansion. 

2.1 Interior Problem: 


The transverse components of fields inside the cavity may be obtained using the 
procedure given in [15]. Expressing the transverse electric and magnetic fields in terms 
of vector modal functions and satisfying the boundary conditions, the fields inside the 
cavity may be written as 


Ei 
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P = 0 


Cos (y (z + L)) » ..(> | -> ^ 

- cfs( T L ~ p ix ’ y) ' f lUi ' c » (x ’ y) r 


2 2 2 2 2 

where y p is the propagation constant, equal to ^Jk - k x - k y , k = k Q e r \i r , ( e r and \i r 

are relative permittivity and permeability of medium inside the cavity ), Y p is the modal 

admittance, e p (x, y) and h p (x, y) are the vector modal functions as defined in [15], 
> 1 — > — > 

j\ = zxH t , H t being the tangential magnetic field over the aperture, and 

-4 —4 

e p (x, y ) = h p (a:, y) xz,z being the unit vector along the z-axis. 

If the 3D surface is divided into triangular subdomains, the electric current over 
the surface of the 3D object including the aperture may be expressed in terms of 
triangular basis functions as [14] 


iV 

){x,y,z ) = ]?T n B n (r') 


th ^ 

where T n is the amplitude of electric current normal to the n edge, B n ( r ') is the 

th 

vector basis fuction associated with n edge, and N is the number of non-boundary edges on the 
surface of the object. The expression for basis function is given by [14] 
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where D n and are the two triangles with the n common edge, A and A are 
+ - — > 

the areas of D n and D w triangles, respectively, and r l and r 2 are position vectors of vertices 

th 4 - 

opposite the the n common edge of and D n triangles, respectively. Likewise the 


magnetic current over the aperture may be written as 
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( 5 ) 


-» M 
m a {x,y,z) = J r m 5 m (r’) 

m = 1 

where T m is amplitude of magnetic current normal to the m edge, and M is the number 
of non-boudary edges over the aperture. Substituting (3) in equations (1) and (2), the total fields 
inside the cavity are obtained as 
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2.2 Electromagnetic Field in Exterior Part: 


( 6 ) 

(7) 


In the exterior part, the total electromagnetic field is obtained by superposing the 

scattered field due to j and m a , and the incident field. 

2.2.1 Incident Field: 


The incident field with time variation e jat 

+ j). 


Ein 


may be written as 

•f * 

-jki • r 


( 8 ) 


where ki - -k Q [xCos ((Jk) Sin (0^) +ySin(ty i )Sin (9 i ) + 2035(0^)] , r = xx + yy + zz, 

Eq = \e\ Cos (a Q ) , and E^ = (a 0 ) , and k Q being the free-space wave number. With 

reference to figure 1, a 0 = 0 corresponds to H-polarization anda 0 = 90° corresponds to E- 
polarization. From equation (8), the x-, y-, and z-components of the incident field may be 
written, respectively, as 

E xi = Eq Cos (0,) Cos (4).) - E^Sin (<j),) (9) 
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Ey- — Eq Cos (0-) Sin ((j)^) + E^Cos 


( 10 ) 


E zi - ~ E e. Sin ( 0 f) 

The corresponding magnetic field components are obtained through 


( 11 ) 


1 > •» 

Hin — T~~ki X Ei n 


k o\ 


( 12 ) 


where r) 0 is the free-space impedance. The incident field with E Q * 0 and E^ = 0 is called 
H-polarized wave and E Q = 0 and E^ * 0 is called the E-polarized wave. 


2.2.2 Scattered Field: 


The scattered field outside the cavity due to ) and m a may be obtained through 
vector electric and magnetic potentials as 



(13) 

m - ivxi 

(14) 

= (-r 0 ) VxF 

(15) 

ro 4e 0 W * f 

(16) 


where the electric and magnetic vector potentials are given by 
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\ds 
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tiYi.^.ff >,£21: 

W 4tc J jU ]>_>, 


-jk 0 {r-r) 


e 0 ,and fi 0 are the permittivity and permeability, respectively, of free-space, and rand 
r' are the coordinates of the field and source points, respectively. Using equations (3) 
and (5), the total scattered electric and magnetic fields may be written as 

E ‘ = l T n- - ^W.A„ J* 2 r m ■ Z ■ < 19 > 


l X~7 

r --VxF m 


N M x 


a)p 0 e 0 


where JV is total number of non-boundary edges over the 3D body including the aperture 
area, and M is total number of non-boundary edges over the aperture area only. A n and F m 

appearing in (19) and (20) are obtained from equations (17) and (18), respectively, by replacing } 

by B n in (17) and m a by B m in (18). 

2.3 Coupled Integral Equations: 

Making the total tangential fields across the aperture continuous and the total 

tangential electric field on the conducting surface of the object to zero we get the following 

coupled integral equations: 

- N / \ M 

2 - 5^ vv4 » J + S • I • vx?„ + = 

Lh = 1 0 m = l 0 J tan 

V -r I X' f J Sin %^ z + L ))\ , N c rf ) ] 

£ n • | X Y ' Cos (y L) e ^ ■ J ‘ e p y ) ) ds j < 21 ) 

n = 1 K P = Q\ P Kl P J J s J 
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N M \ 
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The integral equations given in (21) and (22) may be reduced to a set of algebric equations 

by application of the Method of Moments. By selecting B n ' with n'=\,2,„.N as testing functions 

for equation (21) and B m ' with as testing functions for the equation (22), the 

application of the Method of Moments results in the following matrix equation: 

[c] [z>] [r] = [t] (23) 

>i m lwj m 

The elements of matrices £qj , jjp] , J/0 , and £q| are given by 

= J j B « * ( — jOi'In - J ds 


if ■ (J J • *p <*> y> )* - J J ((X • % <* ))*) 


4 *. = J • I • Vx? "* 


13 



(26) 


Pm'n = J * | ' VxA„ 


ds 
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The elements of column matrices jjyj and [yj are given by 

u n . = ~\\Ein*B n ds (27) 

D* 

V m , = - J \Hin • Bm'ds (28) 

In evaluating numerical values of expressions (24)-(28), numerical integration over 
the triangles are performed using thirteen point gauss quadrature formula [16]. For evaluation of 
self terms; i.e., when n = n' or m = m ' , closed form expressions for these expressions given in 
[17] are used. The unknown electric and magnetic current amplitudes obtained after solving the 
matrix equation (23) can be used to determine scattered far field using expression (19) and (20). 
2.4 Scattered Far Field: 


Using the far field approximation, the scattered electric far field may be obtained from 
equation (19) as 

N M 

(29) 
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E ,e = - X 7 W“ A »e- S T J k o ' F I 
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where A nQ ,A B(j) , F mQ , and F m(J) are given by 


-> e Jk ° r rf jfc 0 Sl 'i 

A„ = 6A„ e -M,A„ + = — 

s„ 


s Sin (9) [^C<?5(4>) + y'Sin ($)] +jk Q z’Cos(d) 


ds n (31) 


^ ^ ~Jk 0 r 

F m = GF mQ + ^ m(fy = 777 j n ”m 

s. 


reft, jk 0 Sin(.B) lx'Cos($) +y'Sin(ty))+jk 0 z'Cos(Q)'\ 

JJ [B m -e )ds n (32) 


The copolarized radar cross section patterns of an aperture backed by a cavity in a 3D conducting 
surface can be obtained from 

J/? J 2 

(33) 
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for H-polarized incidence and 
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for E-polarized incidence. The cross-polarized radar cross section pattern is obtained from 
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3. CODE IMPLEMENTATION 


The coupled field integral equation approach to solve the problem of scattering from a 
rectangular waveguide embedded in a 3D conducting surface has been implemented through a 
code named scatt_reap_recav_3d (scattering from rectangular aperture backed by a rectangular 
casitv embedded in a ^conducting object). 

For running the scatt_reap_recav_3d code, the given geometry is modelled using 
COSMOS/M. In using COSMOS/M, it is assumed that user is familiar with operation of 
COSMOS/M. All dimensions of 3D surface and rectangular cavity are normalized with respect to 
the operating wavelength. 

The input variables must be defined in a file called *.in before running scatt_reap_recav_3cL 
A sample of *. in is, given in Table 1. 


Table 1: Description of* in file 


Variable 

Description 

*.MOD 

Name of input file containing nodes and element infor- 
mation 

output file name 

Name of output file where output is stored 

aa,bb 

x- and y-dimensions of rectangular cavity in wavelength 

alO 

length of cavity in w avelength 

ihigher 

an integer, if zero only dominant mode in the cavity is 
considered. It also skips next line if “ ihigher “ is zero 

ite, itm 

integers, read only if “ ihigher” is not zero, ite is number 
of TE modes and itm is number of TM modes to be con- 
sidered in the cavity 

alpha 

0 for H-polarization 
90 for E-polarization 
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Table 1: Description of* .in file 


Variable 

Description 

theta, phi 

incident theta and phi angles in degrees 

thine, phinc 

increments in theta and phi 

die, emie 

complex relative permittivity and permeability of 
medium inside the cavity 


To demonstrate use of the code, the following examples are considered in this report 

1) rectangular cavity without a ground plane 

2) rectangular cavity embedded in a conducting cube 

3) rectangular cavity embedded at the end of finite circular cylinder 

4) rectangular cavity embedded in the curve surface of finite circular cylinder 


3.1 EXAMPLE 1 

To illustrate use of variables defined in Table I, a rectangular cavity shown in figure 3 is 
taken as an example. To generate *.MOD file for the geometry of figure 3 , fyure3.S‘ES file given 
in Table 1.1 is first run with COSMOS/M. The *.MOD file thus generated is named as 
figure3MOD. The fyure3.in file used for running the scatt_ reap_recav_3d code is shown in Table 1.2 
and a sample output of scatt_reap_recav_3d shown in Table 1.3. 

Table 1.1 Listing of figure3.SES 


C* COSMOS/M Geostar VI. 70 

C* Problem: cube Date: 11-29-94 Time: 8:47:4 

C* 

PT1000 
PT 2 -0.15 -0.15 0 
PT 3 0.15 -0.15 0 
PT 4 0.15 0.15 0 













PT 5 -0.15 0.15 0 
PT 6 -0.15 -0.15 -0.2 
PT 7 0.15 -0.15 -0.2 
PT 8 0.15 0.15 -0.2 
PT 9-0.15 0.15-0.2 
SCALE 0 
CRLINE 1 2 3 
CRLINE 2 3 4 
CRLINE 3 4 5 
CRLINE 4 5 2 
CRLINE 5 6 7 
CRLINE 6 7 8 
CRLINE 7 8 9 
CRLINE 8 9 6 
CRLINE 9 2 6 
CRLINE 10 3 7 
CRLINE 11 4 8 
CRLINE 12 5 9 
CT 1 0 0.06 4 12 3 4 0 
CT 2 0 0.06 4 2 10 6 110 1 
CT 3 0 0.06 4 6 7 8 5 0 1 
CT 4 0 0.06 4 4 9 8 12 0 1 
CT 5 0 0.06 4 1 10 5 9 0 1 
CT 6 0 0.06 4 3 11 7 12 0 1 
RG 1 1 1 0 
RG 2 1 2 0 
RG 3 1 3 0 
RG 414 0 
RG 5 1 5 0 
RG 6 1 6 0 

PH 1 RG 1 0.06 0.0001 1 
MA_PH 111 
PRG 111114 
NMERGE 
NCOMPRESS 
NTCR 1 1 1 1 
NPCR 2221 
NJCR 3 3 3 1 
QCR444 1 

Table 1.2 Listing of fyure3.in 


figure3.MOD “ Input file with node and element information” 
figure3.out “ Output file with RCS as a function of look angle” 

0.3,0.3 “ x- and y- dimensions normalized with wavelength of rectangular cavity 

0.2 length of cavity normalized with respect to wavelength 
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1 integer flag to consider higher order modes in the cavity 

1 0 number of TE modes to be considered 

10 number of TM modes to be considered 

0 alpha = 0 for H-polarization 

0.,0. theta and phi incident angles in degrees 

1 .,0 increment in theta and phi in degrees 

30. frequency in GHz 

(l.,0.),(l.,0.) complex relative permittivity and permeability of medium inside cavity 
Table 1.3 Listing of fyju.re3.out 


X-dimension of WG= 0.3000000 wavelength 

Y-dimension ofWG= 0.3000000 wavelength 

Length of WG cavity = 0.2000000 wavelength 

Number of waveguide modes = 20 

Number of Nodes Used = 128 

Number of Elements Used = 252 

Number of nonboundary edges = 378 

Number of nonboundary edges (aperture) = 65 

Frequency in Ghz = 30.00000 


alphatheta 

phi 

cross-pol(RCS)co-pol(RCS ) 

0 . 0 . 

0 . 

-42.90524 

0.1645633 

0 . 1.000000 

0 . 

- 42.87325 

0.1629779 

0.2.000000 

0 . 

-42.85700 

0.1584223 

0.3.000000 

0 . 

-42.85587 

0.1508932 

0.4.000000 

0 . 

-42.86822 

0.1403691 

0.5.000000 

0 . 

-42.89356 

0.1268443 

0.6.000000 

0 . 

-42.93044 

0.1102918 

0.7.000000 

0 . 

- 42.97783 

9.0677880E-02 

0.8.000000 

0 . 

-43.03438 

6.7987286E-02 

0.9.000000 

0 . 

-43.09773 

4.2179674E-02 

0.10.00000 

0 . 

-43.16571 

1.3209175E-02 

********** 

*** 

******** 

************ 

********** 

*** 

******** 

************ 

0.171.0000 

0 . 

-49.14130 

-0.9871500 

0.172.0000 

0 . 

-49.60229 

-0.9702910 

0.173.0000 

0 . 

-50.03440 

-0.9554753 

0.174.0000 

0 . 

-50.42996 

-0.9426595 

0.175.0000 

0 . 

-50.78228 

-0.9318268 

0.176.0000 

0 . 

-51.08117 

-0.9229562 

0.177.0000 

0 . 

-51.32377 

-0.9160256 

0.178.0000 

0 . 

-51.50285 

-0.9110218 

0.179.0000 

0 . 

-51.61285 

-0.9079247 

0.180.0000 

0 . 

-51.65466 

-0.9067491 



3.2 EXAMPLE 2 


A rectangular cavity embedded in a metal cube as shown in figure 4 is considered as a second 
example. The *.SES, *. in, and *. out files used to for the problem are shown in Tables 2.1, 2.2, 
and 2.3. 


Table 2.1 Listing of figu.re4.SES 


C* 

C* COSMOS/M Geostar VI. 70 

C* Problem : figur4 Date: 1-11-95 Time : 10:25:41 

C* 

PLANE Z 0 1 
VIEW00 1 0 
PT 1 -0.25 -0.25 0 
PT 2 0.25 -0.25 0 
PT 3 0.25 0.25 0 
PT 4 -0.25 0.25 0 
SCALE 0 

PT 5 -0.25 -0.25 -0.5 
PT 6 0.25 -0.25 -0.5 
PT 7 0.25 0.25 -0.5 
PT 8 -0.25 0.25 -0.5 
PT 9 -0.15 -0.15 0 
PT 10 0.15-0.15 0 
PT 110.15 0.15 0 
PT 12-0.15 0.15 0 
VIEW 1110 
SCALE 0 
CRLINE 112 
CRLINE 2 2 3 
CRLINE 3 3 4 
CRLINE 4 4 1 
CRLINE 5 5 6 
CRLINE 6 6 7 
CRLINE 7 7 8 
CRLINE 8 8 5 
CRLINE 9 1 5 
CRLINE 10 2 6 
CRLINE 11 3 7 
CRLINE 12 4 8 
CRLINE 13 9 10 
CRLINE 14 10 11 
CRLINE 15 11 12 
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CRJLINE 16 12 9 
CLS 1 

CT 1 0 0.07 4 1 2 3 4 0 
CT 2 0 0.07 4 13 14 15 16 0 
RG 1 2 1 2 0 
RG 2 1 2 0 

CT 3 0 0.07 4 2 10 6 110 1 
RG 3 1 3 0 

CT 4 0 0.07 4 3 117 12 0 1 
RG 4 14 0 

CT5 0 0.07 4 4 9 8 1201 
RG 5 1 5 0 

CT 6 0 0.07 4 1 10 5 9 01 
RG 6 1 6 0 

CT 7 0 0.07 4 5 6 7 8 0 1 
RG 7 1 7 0 

PH 1 RG 1 0.07 0.0001 1 
MA_PH 1 1 1 
PRG 212114 
NTCR 13 1 13 1 
NPCR 14 2 14 1 
NJCR 15 3 15 1 
QCR 16 4 16 1 


Table 2.2 Listing of fyure4.in 


figure4.MOD 

figure4.out 

0.3,03 

0.2 

1 

30 

30 

0 

0 ., 0 . 

l .,0 

30. 

(l.,0.),(l.,0.) 


“ Input file with node and element information” 

“ Output file with RCS as a function of look angle” 

“ x- and y- dimensions normalized with wavelength of rectangular cavity 

length of cavity normalized with respect to wavelength 

integer flag to consider higher order modes in the cavity 

number of TE modes to be considered 

number of TM modes to be considered 

alpha = 0 for H-polarization 

theta and phi incident angles in degrees 

increment in theta and phi in degrees 

frequency in GHz 

complex relative permittivity and permeability of medium inside cavity 


Table 2.3 Listing of fyure4.out 

X-dimension of WG = 0.3000000 wavelength 
Y-dimension of WG = 0.3000000 wavelength 



Length of WG cavity = 0.2000000 wavelength 
Number of waveguide modes = 60 

Number of Nodes Used = 152 
Number of Elements Used = 300 


Frequency in Ghz = 

30.00000 



alpha 

theta 

phi 

cross-pol 

co-pol 

0 . 

0 . 

0 . 

-13.86164 

3.597274 

0 . 

1.000000 

0 . 

-13.86652 

3.586681 

0 . 

2.000000 

0 . 

-13.87791 

3.557459 

0 . 

3.000000 

0 . 

-13.89583 

3.509618 

0 . 

4.000000 

0 . 

-13.92023 

3.443209 

0. 

5.000000 

0 . 

-13.95116 

3.358291 

0 . 

6.000000 

0 . 

-13.98853 

3.254961 

0 . 

7.000000 

0 , 

-14.03238 

3.133344 

0 . 

8.000000 

0 . 

-14.08267 

2.993612 

0 . 

9.000000 

0 . 

-14.13940 

2.836054 

0 . 

10.00000 

0 . 

-14.20255 

2.660966 

0 . 

11.00000 

0 . 

-14.27210 

2.468850 

0 . 

12.00000 

0 . 

-14.34806 

2.260279 

0 . 

13.00000 

0 . 

-14.43036 

2.036047 

0 . 

14.00000 

0 . 

-14.51905 

1.797124 

0 . 

15.00000 

0 . 

-14.61407 

1.544816 

0 . 

16.00000 

0 . 

-14.71538 

1.280692 

0 . 

17.00000 

0 . 

-14.82301 

1.006683 

0 . 

18.00000 

0 . 

-14.93689 

0.7251824 

0 . 

19.00000 

0 . 

-15.05705 

0.4389912 

0 . 

20.00000 

0 . 

-15.18342 

0.1514122 

0 . 

21.00000 

0 . 

-15.31600 

-0.1337813 

0 . 

22.00000 

0 . 

-15.45476 

-0.4124082 

0 . 

23.00000 

0 . 

-15.59963 

-0.6800060 

0 . 

24.00000 

0 . 

-15.75065 

-0.9319531 

0 . 

25.00000 

0 . 

-15.90774 

-1.163811 






0 . 

170.0000 

0 . 

-32.80526 

2.489367 

0 . 

171.0000 

0 . 

-32.59278 

2.690944 

0 . 

172.0000 

0 . 

-32.40409 

2.870220 

0 . 

173.0000 

0 . 

-32.23733 

3.027478 

0 . 

174.0000 

0 . 

-32.09236 

3.163024 

0 . 

175.0000 

0 . 

-31.97046 

3.277061 

0 . 

176.0000 

0 . 

-31.86998 

3.369819 

0 . 

177.0000 

0 . 

-31.79082 

3.441491 

0 . 

178.0000 

0 . 

-31.73413 

3.492200 

0 . 

179.0000 

0 . 

-31.69735 

3.522068 

0 . 

180.0000 

0 . 

-31.68309 

3.531136 
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3.3 Example 3 

In this example a rectangular aperture backed by a rectangular cavity placed at the one end of a 
finite circular cylinder as shown in figure 5 is considered. *. SES, *. in, and *. out files used for 
this problem are shown in Tables 3.1, 3.2, and 3.3, respectively. 

Table 3.1 Listing of fyure5.SES 


C* 

C* COSMOS/M Geostar VI. 70 

C* Problem : figure5 Date : 1-25-95 Time : 15:52:43 

C* 

C* FILE tempi. SES 1111 
PLANE Z 0 1 
VIEW00 10 
PT100 0 
PT 2 0.35 0.35 0 
SCALE 0 

CRPCIRC 1 1 2 0.494975 360 4 

SCALE 0 

SCALE 0 

PT 6 0 0 -0.5 

PT 7 0.35 0.35 -0.5 

SCALE 0 

VIEW 1110 

CRPCIRC 5 6 7 0.494975 360 4 

PT 11 -0.15-0.15 0 

PT 12 0.15-0.15 0 

PT 13 0.15 0.15 0 

PT 14-0.15 0.15 0 

CRLINE 9 2 7 

CRLINE 10 3 8 

CRLINE 114 9 

CRLINE 12 5 10 

CRLINE 13 11 12 

CRLINE 14 12 13 

CRLINE 15 13 14 

CRLINE 16 14 11 

CRLINE 17 12 5 

CRLINE 18 13 2 

CRLINE 19 14 3 

CRLINE 20 1 1 4 
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SF4CR 1 3 17 13 20 0 
SF4CR 2 17 4 18 14 0 
SF4CR 3 18 1 19 15 0 
SF4CR4 19 2 20 16 0 
SF4CR 5 13 14 15 16 0 
SF4CR 6 1 9 5 10 0 
SF4CR7 10 6112 0 
SF4CR8 113 12 7 0 
SF4CR 9 12 8 9 4 0 
SF4CR 10 5 6 7 8 0 
PH 1 SF 1 0.1 0.0001 1 
MA_PH 1 1 1 
NTCR 13 1 13 1 
NPCR 14 2 14 1 
NJCR 15 3 15 1 
QCR 16 4 161 
PSF 5 1 5 1 1 1 4 
NMERGE 1 831 1 0.0001 0 1 0 
NCOMPRESS 1 831 


Table 3.2 Listing oifyureS.in 


figure5.MOD 

figure5.out 

0.3,03 

0.2 

1 

10 

10 

0 

0 ., 0 . 

L,0 

30. 


“ Input file with node and element information” 

“ Output file with RCS as a function of look angle” 

“ x- and y- dimensions normalized with wavelength of rectangular cavity 

length of cavity normalized with respect to wavelength 

integer flag to consider higher order modes in the cavity 

number of TE modes to be considered 

number of TM modes to be considered 

alpha = 0 for H-polarization ( for E-polarization alpha =90) 

theta and phi incident angles in degrees 

increment in theta and phi in degrees 

frequency in GHz 

complex relative permittivity and permeability of medium inside cavity 


Table 3.3 Listing of fyure5. out 

X-dimensionofWG = 0.3000000 wavelength 
Y-dimensionofWG= 0.3000000 wavelength 
Length of WG cavity = 0.2000000 wavelength 
Number of waveguide modes = 20 

Number of Nodes Used = 831 
Number of Elements Used = 1658 

Frequency in Ghz = 30.00000 
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alpha 

theta 

phi 

cross-pol 

co-pol 

0 . 

0 . 

0 . 

-5.916269 

9.834730 

0 . 

1.000000 

0 . 

-5.920516 

9.820421 

0 . 

2.000000 

0 . 

-5.933593 

9.778456 

0 . 

3.000000 

0 . 

-5.955378 

9.709249 

0 . 

4.000000 

0 . 

-5.985921 

9.613455 

0 . 

5.000000 

0 . 

-6.025254 

9.492058 

0 . 

6.000000 

0 . 

-6.073274 

9.346359 

0 . 

7.000000 

0 . 

-6.129971 

9.177999 

0 . 

8.000000 

0 . 

-6.195318 

8.988997 

0 . 

9.000000 

0 . 

-6.269305 

8.781787 

0 . 

10.00000 

0 . 

-6.351826 

8.559145 

0 . 

11.00000 

0 . 

-6.442932 

8.324296 

0 . 

12.00000 

0 . 

-6.542520 

8.080740 

0 . 

13.00000 

0 . 

-6.650455 

7.832253 


0 . 

169.0000 

0 . 

-33.84848 

7.762480 

0 . 

170.0000 

0 . 

-33.33032 

8.084915 

0 . 

171.0000 

0 . 

-32.86591 

8.382092 

0 . 

172.0000 

0 . 

-32.45417 

8.652040 

0 . 

173.0000 

0 . 

-32.09342 

8.893197 

0 . 

174.0000 

0 . 

-31.78228 

9.104296 

0 . 

175.0000 

0 . 

-31.52011 

9.284355 

0 . 

176.0000 

0 . 

-31.30555 

9.432592 

0 . 

177.0000 

0 . 

-31.13779 

9.548450 

0 . 

178.0000 

0 . 

-31.01659 

9.631517 

0 . 

179.0000 

0 . 

-30.94126 

9.681516 

0 . 

180.0000 

0 . 

-30.91166 

9,698285 


3.4 Example 4 

A rectangular aperture backed by a rectangular cavity and placed over the curved surface of 
finite circular cylinder as shown in figure 6 is considered. *. SES, *. in, and *. out files used 
this case are shown in Tables 4.1, 4.2, and 4.2, respectively. 


Table 4.1 Listing of fyure6.SES 


C* 

C* COSMOS/M Geostar VI. 70 

C* Problem : figure6 Date: 1-26-95 Time: 8:10:19 



c* 

C* FILE tempi. SES 1111 
PLANE Z 0 1 
VIEW 0 0 10 
PT 1 0 0 0 

PT 2 -0.35 -0.1505 0 
PT 3 0.35 -0.1505 0 
PT 4 0.35 0.1505 0 
PT 5 -0.35 0.1505 0 
PT 6 -0.35 -0.1505 0.1 
PT 7 0.35 -0.1505 0.1 
PT 8 0.35 0.1505 0.1 
PT 9 -0.35 0.1505 0.1 
SCALE 0 
VIEW 1110 
PT 10 -0.35 0 -0.38 
SCALE 0 
PLANE X 0 1 

CRPCIRC 1 10 9 0.503041 34.82 1 

CRPCIRC 2 10 6 0.503041 55.18 1 

CRPCIRC 3 10 11 0.503041 270 3 

SCALE 0 

PT 14 -0.5 0 -0.38 

PT 15-0.5 0.1505 0.1 

PT 16-0.5 -0.1505 0.1 

CRPCIRC 6 14 15 0.503041 34.82 1 

CRPCIRC 7 14 16 0.503041 55.18 1 

CRPCIRC 8 14 17 0.503041 270 3 

CT 100.15 5 67 8 9 100 

RG 1 1 1 0 

CRLINE 11 15 9 

CRLINE 12 16 6 

CRLINE 13 17 11 

CRLINE 14 18 12 

CRLINE 15 19 13 

SF4CR 1 6 12 1 11 0 

SF4CR 2 12 7 13 2 0 

SF4CR 3 13 3 14 8 0 

SF4CR 4 14 4 15 9 0 

SF4CR 5 15 10 11 5 0 

CRLINE 16 9 5 

CRLINE 17 5 2 

CRLINE 18 6 2 

SF4CR 6 1 18 17 16 0 

CRLINE 19 9 8 

CRLINE 20 5 4 
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CRLINE 21 8 4 
SF4CR 7 16 19 21 20 0 
CRLINE 22 8 7 
CRLINE 23 4 3 
CRLINE 24 7 3 
PT 20 0.35 0-0.38 
PLANE X 0 1 

CRPCIRC 27 20 8 0.503041 34.82 1 

SF4CR 8 27 24 23 21 0 

CRLINE 25 6 7 

CRLINE 26 2 3 

SF4CR 9 18 26 24 25 0 

SF4CR 10 17 20 23 26 0 

CRPCIRC 28 20 7 0.503041 55.18 1 

CRPCIRC 29 20 21 0.503041 270 3 

CRLINE 32 9 6 

PT 24 0.5 0 -0.38 

PT 25 0.5 0.1505 0.1 

PT 26 0.5-0.1505 0.1 

PLANE X0 1 

CRPCIRC 37 24 25 0.503041 34.82 1 

CRPCIRC 38 24 26 0.503041 55.18 1 

CRPCIRC 39 24 27 0.503041 270 3 

CT 2 0 0.15 5 37 38 39 40 410 

RG 2 1 2 0 

CRLINE 42 8 25 

CRLINE 43 7 26 

SF4CR 13 42 37 43 27 0 

CRLINE 44 21 27 

SF4CR 14 43 38 44 28 0 

CRLINE 45 22 28 

SF4CR 15 29 44 39 45 0 

CRLINE 46 23 29 

SF4CR 16 45 40 46 30 0 

SF4CR 17 46 41 42 31 0 

CRLINE 47 1121 

CLS 1 

SF4CR 18 2 47 28 25 0 
CRLINE 48 12 22 
SF4CR 19 3 47 29 48 0 
CRLINE 49 13 23 
SF4CR 20 4 48 30 49 0 
SF4CR 21 5 49 31 19 0 
PH 1 SF 1 0.1 0.0001 1 
CLS 1 

MA_PH 1 1 1 
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/ 


VIEW00 10 
VIEW 1110 
PSF 10 1 10 1 1 1 4 
CLS 1 
CLS 1 

NTCR 26 1 26 1 
NPCR 23 2 23 1 
NJCR 20 3 20 1 
QCR 17 4 17 1 
CLS 1 


Table 4.2 Listing of jtgureS.in 


figure6.MOD 

figure6.out 

0.7,031 

0.1 

1 

30 

30 

0 

0 ., 0 . 

l .,0 

30. 

(l.,0.),(l.,0.) 


“ Input file with node and element information” 

“ Output file with RCS as a function of look angle” 

“ x- and y- dimensions normalized with wavelength of rectangular cavity 

length of cavity normalized with respect to wavelength 

integer flag to consider higher order modes in the cavity 

number of TE modes to be considered 

number of TM modes to be considered 

alpha = 0 for H-polarization ( for E-polarization alpha =90) 

theta and phi incident angles in degrees 

increment in theta and phi in degrees 

frequency in GHz 

complex relative permittivity and permeability of medium inside cavity 


Table 4.3 Listing of figweS.out 

X-dimension of WG = 0.7000000 wavelength 
Y-dimension of WG = 0.3100000 wavelength 
Length of WG cavity = 0.1000000 wavelength 
Number of waveguide modes = 60 

Number of Nodes Used = 552 
Number of Elements Used = 1100 

Frequency in Ghz = 30.00000 


alpha 

theta 

phi 

cross-pol 

co-pol 

0 . 

0 . 

0 . 

-5.933923 

6.347271 

0 . 

1.000000 

0 . 

-5.935130 

6.339834 

0 . 

2.000000 

0 . 

-5.946536 

6.312681 
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0 . 

3.000000 

0 . 

-5.968001 

6.266007 

0 . 

4.000000 

0 . 

-5.999329 

6.200232 

0 . 

5.000000 

0 . 

-6.040510 

6.115752 

0 . 

6.000000 

0 . 

-6.091329 

6.013002 

0 . 

7.000000 

0 . 

-6.151697 

5.892314 

0 . 

8.000000 

0 . 

-6.221767 

5.753819 

0 . 

9.000000 

0 . 

-6.301365 

5.597322 

0 . 

10.00000 

0 . 

-6.390241 

5.422188 

0 . 

174.0000 

0 . 

-58.31413 

4.871758 

0 . 

175.0000 

0 . 

- 57.79329 

5.028406 

0 . 

176.0000 

0 . 

- 57.26915 

5.159277 

0 . 

177.0000 

0 . 

-56.59153 

5.262630 

0 . 

178.0000 

0 . 

-55.61058 

5.337103 

0 . 

179.0000 

0 . 

-54.31091 

5.381660 

0. 

180.0000 

0 . 

-52.79220 

5.395699 


3.5 Example 5 

In this example a rectangular aperture by backed by a rectangular cavity and placed along the 
circumference of a finite conducting circular cylinder, as shown in figure 7 is considered. *.SES, 
*. in, and *. out files used to run this case are shown in Tables 5.1, 5.2, and 5.3, respectively. 


Table 5.1 Listing of fyure7.S‘ES 
C* 

C* COSMOS/M Geostar V1.70 

C* Problem : figure7 Date : 1-30-95 Time : 11:33:48 

C* 

C* FILE figure7.SES 1111 

PLANE Z 0 1 

VIEW 0 0 10 

PT 1 0 0 0 

PT 2 -0.1505 -0.35 0 

PT 3 0.1505 -0.35 0 

PT 4 0.1505 0.35 0 

PT 5 -0.1505 0.35 0 
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PT 6 -0.1505 -0.35 0.1 
PT 7 0.1505 -0.35 0.1 
PT 8 0.1505 0.35 0.1 
PT 9 -0.1505 0.35 0.1 
PT 10 -0.1505 0 -0.25 
SCALE 0 
VIEW 1110 
PLANE X 0 1 

CRPCIRC 1 10 9 0.494975 90 1 
CRPCIRC 2 10 6 0.494975 270 3 
SCALE 0 

PT 13-0.1505 0-0.25 

PT 13 -0.5 0 -0.25 

PT 14 -0.5 0.35 0.1 

PT 15 -0.5 -0.35 0.1 

CRPCIRC 5 13 14 0.494975 90 1 

CRPCIRC 6 13 15 0.494975 270 3 

SCALE 0 

CT 100.1 4 5 67 8 0 

RG 1 1 1 0 

CRLINE 10 15 6 

CRLINE 11 14 9 

SF4CR 1 1 115 10 0 

CRLINE 12 16 11 

SF4CR 2 10 2 12 6 0 

CRLINE 13 17 12 

SF4CR 3 12 3 13 7 0 

SF4CR4 13 4 11 8 0 

CRLINE 14 9 5 

CRLINE 15 6 2 

CRLINE 16 2 5 

CRLINE 17 5 4 

CRLINE 18 4 3 

CRLINE 19 3 2 

CRLINE 20 8 4 

CRLINE 21 7 3 

PT 18 0.1505 0-0.25 

CRPCIRC 22 18 8 0.494975 90 1 

CRPCIRC 23 18 7 0.494975 270 3 

CRLINE 26 9 8 

CRLINE 27 12 20 

SF4CR 5 26 25 27 4 0 

CRLINE 28 11 19 

SF4CR 6 27 24 28 3 0 

CRLINE 29 6 7 

SF4CR 7 28 23 29 2 0 
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SF4CR91 141615 0 
SF4CR 10 26 20 17 14 0 
SF4CR 11 22 21 18 20 0 
SF4CR 12 29 21 19 15 0 
SF4CR 13 16 19 18 17 0 
PT 21 0.5 0 -0.25 
PT 22 0.5 0.35 0.1 
PT 23 0.5 -0.35 0.1 
SCALE 0 

CRPCIRC 30 21 22 0.494975 90 1 

CRPCIRC 31 21 23 0.494975 270 3 

CT 2 0 0.1 4 30 31 32 33 0 

RG 2 1 2 0 

CRLINE 34 8 22 

CRLINE 35 7 23 

SF4CR 14 34 30 35 22 0 

CRLINE 36 19 24 

SF4CR 15 35 31 36 23 0 

CRLINE 37 20 25 

SF4CR 16 36 32 37 24 0 

SF4CR 17 37 33 34 25 0 

PH 1 SF 1 0.1 0.0001 1 

CLS 1 

MA_PH 1 1 1 
CLS 1 

PSF 13 1 13 1 1 1 4 
NTCR 19 1 19 1 
NPCR 18 2 18 1 
NJCR 17 3 17 1 
QCR 164 16 1 

Table 5.2 Listing of figure7.in 


figure7.MOD 

figure7.out 

0.31,0.7 

0.1 

1 

30 

30 

0 

0 ., 0 . 

l.,0 

30. 

(l.,0.),(l.,0.) 


“ Input file with node and element information” 

“ Output file with RCS as a function of look angle” 

“ x- and y- dimensions normalized with wavelength of rectangular cavity 

length of cavity normalized with respect to wavelength 

integer flag to consider higher order modes in the cavity 

number of TE modes to be considered 

number of TM modes to be considered 

alpha = 0 for H-polarization ( for E-polarization alpha =90) 

theta and phi incident angles in degrees 

increment in theta and phi in degrees 

frequency in GHz 

complex relative permittivity and permeability of medium inside cavity 



Table 5.3 Listing of jyu.Te7.out 

X-dimension ofWG= 0.31000000 wavelength 
Y-dimensionofWG= 0.700000 wavelength 
Length of WG cavity = 0.1000000 wavelength 
Number of waveguide modes = 60 

Number of Nodes Used = 550 
Number of Elements Used = 1 096 

Frequency in Ghz = 30.00000 


alpha 

theta 

phi 

cross-pol 

co-pol 

0 . 

0 . 

0 . 

-14.81105 

6.143041 

0 . 

1.000000 

0 . 

-14.86441 

6.125021 

0 . 

2.000000 

0 . 

-14.92570 

6.073569 

0 . 

3.000000 

0 . 

-14.99266 

5.989009 

0 . 

4.000000 

0 . 

-15.06429 

5.871937 

0 . 

5.000000 

0 . 

-15.13739 

5.723056 

0 . 

6.000000 

0 . 

-15.21013 

5.543458 

0 . 

7.000000 

0 . 

-15.28020 

5.334379 

0 . 

8.000000 

0 . 

-15.34564 

5.097300 

0 . 

9.000000 

0 . 

-15.40529 

4.833958 

0 . 

10.00000 

0 . 

-15.45568 

4.546342 

0 . 

11.00000 

0 . 

-15.49823 

4.236658 

0 . 

12.00000 

0 . 

-15.53147 

3.907334 

0 . 

13.00000 

0 . 

-15.55477 

3.560932 






0 . 

172.0000 

0 . 

-35.59557 

4.241209 

0 . 

173.0000 

0 . 

-35.00657 

4.387676 

0 . 

174.0000 

0 . 

-34.46977 

4.516529 

0 . 

175.0000 

0 . 

-33.97979 

4.627316 

0 . 

176.0000 

0 . 

-33.53917 

4.719339 

0 . 

177.0000 

0 . 

-33.14160 

4.791856 

0 . 

178.0000 

0 . 

-32.78939 

4.844146 

0 . 

179.0000 

0 . 

-32.47799 

4.875624 

0 . 

180.0000 

0 . 

-32.20755 

4.885905 
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4. NUMERICAL RESULTS 

To validate the present analysis and computer code, various numerical examples are 
considered. For validation, the numerical results obtained by ME-MoM are compared with 
results obtained by the hybrid Finite Element Method and Method of Moments ( FEM-MoM ) 
method. In the FEM-MoM method, EM fields inside the cavity are obtained using the Finite 
Element Method and the fields outside the cavity are obtained using the Method of Moments. 
The unknown fields are obtained by subjecting the total electric fields to appropriate boundary 
conditions. The results obtained by ME-MoM are also compared with the pure MoM approach. 
In the pure MoM approach, the inner wall of cavity along with the outer surface of the conducting 
body is considered as one closed 3D conducting surface. The Electric Field Integral Equation 
(EFEE) with electric surface current density as unknown is then obtained by subjecting the total 
tangential electric field on the entire closed surface to zero. The electric surface current density 
obtained after solving the EF1E is then used to determine scattering pattern of cavity backed 
aperture recessed in a 3D conducting body. 

As a first example, the radar cross section of a rectangular cavity, as shown in figure 
3, with dimensions a= 0.3 X 0 , b=0.3A. 0 and c = 0.2X 0 is calculated. The cavity is assumed to be 
open at the z = 0 plane. For the numerical solution of equation (23), the infinite summations with 
respect to index p must be truncated to some finite number P , where P is the number of modes 
considered in ascending order of their cutoff frequencies. Numerical convergence of matrix 
equation (23) depends upon the choices of values of M, IV, and P . For numerical convergence, 
sufficiently large values of M, N, and P must be selected. For M = 42 , N = 252, and 
P = 20 , the Radar Cross Section (RCS) of the rectangular cavity is calculated and presented in 
figure 3 along with the results obtained from the hybrid FEM-MoM. Further increases in values 


33 



M , N , and P were found to result in very insignificant changes in the RCS. The numerical 
results obtained using both methods agree very well verifying the validity of the combined field 
integral equation approach presented in this paper. 


A second example considered for validation of the present approach is a rectangular 
cavity embedded in a solid conducting cubical box as shown in figure 4. The cavity dimensions 
are the same as described in figure 3. The cubical box with sides equal to 0.5 A is considered. 
For N = 300 , M = 42 , and P = 60 , the RCS pattern of the cavity embedded in a cubical solid 
box is calculated and presented in figure 4 along with the numerical results obtained using the 
pure MoM approach. In the pure MoM approach, the cubical box with the inner surface of 
cavity was considered as a closed conducting surface. From the plots in figure 4, it may be 
concluded that the results obtained by both methods agree well. 

The present approach is also used to predict monostatic RCS’s of cavities embedded 
in a finite metallic circular cylinder. In figure 5, a monostatic pattern of a rectangular cavity with 
dimensions a = 0.3 X, b = 0.3 A, and c == 0.2 A,, which is embedded in one of the ends of the 
finite circular cylinder, is presented. The finite metallic cylinder is of length L - XX and 
diameter r - 0.5 X. The monostatic pattern of the cavity is also calculated using the pure MoM 
approach and is presented in figure 5. The results obtained by both methods agree well. 

To demonstrate the application of the present approach to a cavity backed aperture 
mounted on curve surfaces, a rectangular cavity withx width equal to 0.7 X, y width equal to 
0.31 A, and z depth equal to 0.2 A embedded in the finite metallic cylinder as shown in figure 6 is 
considered. To be able to match the cavity modal fields with electromagnetic fields outside the 
cavity across the cavity aperture, the aperture surface must be planar and normal to z-axis. 
However, for the example shown in figure 6 the aperture plane is not planar. To avoid this 
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difficulty for apertures on curve surfaces, the equivalent cavity aperture which is a little inside 
the curved surface is considered for matching fields across the cavity aperture. For the example 
under consideration, the cavity aperture at z - -0.1 A is considered. The monostatic pattern of 
the cavity mounted on cylindrical surface calculated using the present method and the pure MoM 
method is shown in figure 6. There is good agreement between the two results. 

Monostatic patterns of a rectangular cavity with longer dimension along the 

circumferential direction of a finite metallic cylinder are calculated using the present method 
and the pure MoM technique. The results are presented in figure 7. For these results, a 
rectangular cavity with x width equal to 0.31 A, , y width equal to 0.7 A, , and z depth equal to 
0.2 A. embedded in a finite metallic cylinder of length L equal to 1 .0 A and the diameter r equal to 
0.5 A was considered. The results obtained by both methods agree well with little disagreement 
near broad side angles. 

5. CONCLUSION 

The problem of electromagnetic scattering from a cavity backed aperture recessed in a 
3D conducting body has been analyzed using the combined field integral equation approach. The 
EM fields outside the cavity are obtained in terms of the free space Green’s function and electric 
and magnetic surface current densities present on the surface of conducting bodies including the 
aperture. The EM fields inside the cavity are obatined in terms of cavity modal functions. The 
combined field integral equations are then derived by subjecting the electric and magnetic fields 
to appropriate boundary conditions. Using the Method of Moments, the integral equations are 
solved for unknown surface current densities. The scattering characteristic of the cavity backed 
aperture is then determined from the surface current densities. The numerical results obtained by 
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the present method are compared with the numerical results obtained by FEM-MoM and the pure 
MoM methods. The computed results for the rectangualr cavities considered in this paper agree 
well with the results obtained by the hybrid FEM-MoM and the pure MoM methods. The 
limitation of the present method is that it can be used to analyze regularly shaped cavities filled 
with homogeneous material. 
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Backscatter RCS patterns for a rectangular cavity as shown with dimensions 
a=0.3 X , b= 0.3 X , c= 0.2 X without a ground plane for E- and H- polarized 
incident plane wave. Solid and hollow triangles indicate numerical data 
obtained using FEM-MoM 
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Backscatter RCS patterns for a rectangular cavity (a=0.3 X , b=0.3 X, c=0.2X ) 
embedded in a solid cube with side S =0.5 X , as shown 

for E- and H-polarized plane wave incidenc ( <J>j = 0°). Solid and hollow triangles 
indicate results obtained by the method of moments. 
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Figure 6 Backscatter RCS pattern of rectangular cavity with x width = 0.7 X, y width 
0.31 X, and depth = 0.2 X, embedded in a conducting circular cylinder with L 
= 1 X. r = 0.5 X. 
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Figure 7 Backscatter RCS pattern of rectangular cavity with x width = 0.31 A, y 

width =0.7X, and depth = 0.2 X embedded in a conducting circular cylind 
with L = 1 .OX. and r = 0.5X. 
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